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ABSTRACT 

We perform a set of non-radiative hydrodynamical simulations of merging spherical halos in order to un- 
derstand the angular momentum (AM) properties of the galactic halos seen in cosmological simulations. The 
universal shape of AM distributions seen in simulations is found to be generically produced as a result of 
mergers. The universal shape is such that it has an excess of low AM material and hence cannot explain the 
exponential structure of disk galaxies. A resolution to this is suggested by the spatial distribution of low AM 
material which is found to be in the center and a conical region close to the axis of rotation. A mechanism 
that preferentially discards the material in the center and prevents the material along the poles from falling 
onto the disc is proposed as a solution. We implement a simple geometric criteria for selective removal of low 
AM material and show that in order for 90% of halos to host exponential discs one has to reject at least 40% 
of material. Next, we explore the physical mechanisms responsible for distributing the AM within the halo 
during a merger. For dark matter there is an inside-out transfer of AM, whereas for gas there is an outside-in 
transfer, which is due to differences between collisionless and gas dynamics. This is responsible for the spin 
parameter A and the shape parameter a of AM distributions being higher for gas as compared to dark matter. 
We also explain the apparent high spin of dark matter halos undergoing mergers and show that a criteria stricter 
than what is currently used, would be required to detect such unrelaxed halos. Finally, we demonstrate that the 
misalignment of AM between gas and dark matter only occurs when the intrinsic spins of the merging halos are 
not aligned with the orbital AM of the system. The self-misalignment (orientation of AM when measured in 
radial shells not being constant), which could be the cause of warps and anomalous rotation in disks galaxies, 
also occurs under similar conditions. The frequency and amplitude of this misalignment is roughly consistent 
with the properties of warps seen in disk galaxies. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: halos — galaxies: kinematics and dy- 
namics 



1. INTRODUCTION 

In the Standard picture of galaxy formation, galactic ha- 
los acquir e their angular momentum (hereafter AM) via tidal 
torques ( Peeblesl [r969) in the linear regime and the process 
lasts till about turnaround, when the system decouples from 
the Hubble flow. After the collapse the system forms a virial- 
ized Structure. The gas inside the virialized dark matter halo 
then cools radiatively and collapses while conserving its AM, 
resulting in the formation of centrifugally supported disks 
(IWhite & ReesI [T978I: iFafl & Efstathlpul [T98QI: fWhita 1984). 
The process is also accompa nied by the adiabatic c ontrac- 
tion of the dark matter halo (iBlumenthal et al.iri986h . This 
Standard picture leads to distribution of size and luminos- 
ity of galaxies in reas onable agreement with observations 
JKauffmann 1996; Da lcanton etaP 119971: IMo et alJ [T99l; 
lAvila-Reese et al. 1998;; Ivan den BoschI I200QI: iDutton et al.1 
|2007; Gnedin et al. 2007]). 

But detailed simulations revealed two problems. Firstly 
in simulations incorporating gas with cooling and star 
formation, the gas was found to lose a significant fraction 
of its AM, resulting in disks which were too small in 
size, a problem known a s the a ngular momentun i catas- 
trophe ("Navarro & Benz' '1991^, 'Navarro & White' '1994^, 
[Navarro & Steinmetz 1997; Steinmetz & Navarro 1999; 
ISommer-Larsen et al. 1999). The cause of the problem being 
that due to efficient cooling, the gas is accreted as dense 
clumps which during mergers loses its AM via dynamical 



friction. 

Second problem is the angular momentum distribution 
(AMD hereafter) problem, i.e., even if the AM is as- 
sumed to be conserved one cannot explain the exponen- 
tial nature of disk galaxies. Using cold dark mat t er nu- 
merical simulations, it was shown by iBullock et al.l (1200 Ih 
that if disks are formed from gas with AMDs similar to 
that of dark matter, then this results in excess mass near 
the center as compared to an exponential disc. Specifi- 
cally there is too much low AM material and this makes 
it v ery hard to explain the ori g in of bulgeless dwar f galax- 
ies ([vandinBoscLilSI ISSoIl; |va^ Sim- 
ulations inc orporation no n-radiative gas also lead to sim ilar 
concl usions (Ivan den Bos ch et al. 2002; Sharma & Steinniet3 
12005^). As demonstrated in Sharma & Steinmetz (2005), the 
resulting AMDs written in terms of s = j/jtot closely fol- 
lows a law of form P{s) = [a^ /r{a)]s^~^e~^^ , the uni- 
versal form found in dark matter halos of cosmological A^- 
body simulations. Although the a parameter for gas is slightly 
higher (close to 0.9) than that of dark matter (0.83), it is still 
much less than what is needed for explaining the exponential 
structure of galactic disks, i.e. , a > 1.3. 

The origin of the the universal form of the AMDs is still 
poorly understood and if we can under stand it, that m a y pro- 
vide the clue to solving the problem. Mailer & Dekell (120021) 
proposed a model of build up of AM by a sequence of merg- 
ers. In this model, the final halo spin is assumed to be the sum 



of orbital angular momenta of merging satellites. The model 
was found to correctly repr oduce the distribution of spi n pa- 
rameters of halos ( Vitvitska erani2002tlMaller et al.ll200 2). A 
simple extension of this model was also found to roughly re- 
produce the AMDs. According to this model, the magnitude 
and direction of the total AM of a halo is predominantly de- 
termined by the last major merger and hence the major merger 
contributes to the high AM part of the AMD. The numerous 
small satellites fall in from random directions and mainly con- 
tribute to the low AM part of the AMD. This suggests that 
blowout of gas, e.g., by means of supernova feedback, from 
small halos can eliminate the low AM part of the distribution 
and may resolve the AMD problem in addition to the angular 
momentum catastrophe. 

An alternative solution to the AMD problem is that the 
feedback driven outflows preferentially disca rd low AM 
materi al during the assembly of the galaxy (|Brook et al.l 
1201 Ibl). In fact recent high resolution simulations including 
star formation and feedback have be en quite successful in 
forming bulge-less exponential disks (iGovemato et al.l 120101: 



iBrooket aLOOllbc iGuedes etal. 1120 lib where such a process 
has been shown to occur. Understanding the spatial distribu- 
tion of the low AM material may tell us as to which method 
is more effective in solving the AMD problem.^ 

According to the model proposed by Mailer & Dekell 
([2002), the most favorable scenario for galaxy formation is, 
where there are very few minor mergers, e.g., a halo acquir- 
ing its AM via a major merger. Is it enough to generate AMDs 
such that exponential disks can be formed? If the gas distribu- 
tion is concentrated due to cooling or puffed up as with feed- 
back, does it change the AMD of merger remnants? These are 
some of the questions that we investigate. 

The AM properties of galaxies is of increasing inter- 
est in observational surveys. For example, the Calar Alto 
Legacy Integral Field Area survey (CALIFA) will obtain spa- 
tially resolved spectroscopic information for about 600 galax- 
ies in the local univ erse, using integral field spectroscopy 
(ISanchez et al.ll20lTh . New imaging fiber bundles (so called 
hexabundle s) are to be used on wide-fie l d survey t elescopes 
(e.g., AAT: iBland-Hawthorn et ani201 ll: IcFoom e t al. 2011) 
to obtain spatially resolved stellar and gas kinematics for a 
volume-limited sample 10^~^ galaxies. With such surveys it 
will be possible to study the AMD of galaxies in voids, fil- 
aments, groups and clusters. Although current simulations 
which include star formation and feedback have started show- 
ing success in forming disk galaxies, but these simulations 
are computationally very expensive and this prohibits gen- 
eration of a large sample of galaxies for statistical studies. 
On the other hand, dark matter only simulations are com- 
putationally much less demanding which makes them suit- 
able for comparison with large scale galaxy surveys, but one 
needs a way to populate dark matter halos with galaxies. 
Semi-analytic modelling of galaxies p rovides a way to do this 



is well fit by a log normal distribution (Ivan den B osch et aL 



(lColeetal .1119941: iBaugh et al.lll996l; feauffmann et^aD 



Somerville & Primack 1999; Kauffmann & Haehnelfl 



1999; 



120001: 



Benson et al.ii2003i: iKang et alj i2005). However, most semi- 



analytic models do not take the low AM problem into account 
and this can have important consequences. Additionally, it is 
often assumed that the AM properties of gas is same as that of 
dark matter. This provides a motivation for studying the dif- 
ferences between the AM properties of gas and dark matter. 

The two main AM properties are the spin parameter and the 
shape of AMDs. The spin of dark matter halos has been ex- 
tensively studied and it has been shown that the distribution 



2002; Sharma & Steinmetz 2005; Bett et al. 2007; N eto et al. 
i2007i: iMaccio et al.i i2007i: iBett et al.i i201Q) . In compari- 
son th ere have been far fewer stud ies on AM properties of 
gas. In Sharma & Steinmetz (2005) and lChen et aH ([2003), it 
was found that for halos simulated in cosmological context, 
the spin parameter and shape parameter of AMDs is higher 
for gas as compared to tha t of dark matter. Additionally, 
ISharma & Steinmet3 (l2005h found the spin ratio Acas/ADM 
to increase with cosmic time (at redshift zero the value be- 
ing 1.4). The cause of this trend is still not known. In 
contrast, van den Bosch etal] (1200 2) found the AM proper- 
ties of gas and dark matter to be very similar. A number of 
reason could be res ponsible for this di screpancy. Firstly, it 
could be because ivan den Bosch et aO (i2002j) in their anal- 
ysis had included a la rge number of halos with low parti- 
cle numbers. Secondly, van d en Bosch et al . (2002) had used 
thermally broadened gas velocities to compare the AMDs of 
gas with that of dark matter, whereas broadening of veloc- 
ities is known to ma s k out the differe nces between AMDs 
ISharma & Steinmet3 (l2005h . Finally, Ivan den Bosch et al.l 
(2002) had analyzed the results at z = 3 whereas the other 
authors had analyzed them at z = 0. 

Recently, it has been reported that high spin halos 
are more clustered than low spin halos (iBett et al.l I2007L 
iDavis & Natarajanll201Qh . iMaccio et aP (120071) on the other 
hand do not find any environmental dependence. A cru- 
cial difference in the two schemes is the treatment of unre- 
laxed halos. It has been shown that out-of-equilibrium ha- 
los tend to have higher spin and lower concentration, which 
whe n removed makes the halo concentration independent of 
spin ('GardneiJ'2001^, 'Vitvitska et al."2002^, 'P eirani et al.ll2004l: 
Hetznecker & Burkert 2006; Neto et al. 20071). Such an effect 
could also be responsible for higher clustering of high spin ha- 
los. D'Onghia & Navarro (2007) have studied the correlation 
of merger history and spin of halos and found that halos im- 
mediately after merging have higher spin. Later on during the 
virialization process the halos spin down due to redistribution 
of mass and AM. Generally the offset of the center of mass 
is used to parameterize the unrelaxed halos. How effective 
is this parameter in detecting unrelaxed halos? Observation- 
ally it is the spin of the baryonic component that is observed, 
hence it is important to know if the gas also undergoes such a 
spin up and spin down during mergers? 

Another area where gas shows a difference from dark mat- 
ter is the issue of misalignment between them. In non radia- 
tive hydrodynamical simulations the AM of gas in galactic 
halos is found to be misaligned with respect to dark mat- 
ter with a mean angle of 20° (Ivan den Bosch et aT] 120021: 
ISharma & Steinmet3l2005|) . In simulations with star forma- 
tion and feedback the galactic discs are also found to be mis- 
ahgned, with a median angle of - 30° (Bett et aLllMTol) . The 
misalignment has important observational consequences. For 
example it has been found that the distribution of satellite 
galaxies is preferentially align ed along the major axis of the 
central galaxy (Brainerd 2005; Yang et al."2006, 'Az zaro et"an 
2007; Wang et al. 2008). Agustsson & Brainerd (2006|) show 
that if the disk AM vectors are aligned with the minor axis of 
the halo or the AM o f the halo th en the observed anisotropy 
can be reproduced. iKang et all (i2007) further showed the 
second option is preferred as orientation with minor axis re- 
sults in a stronger signal than that observed. If the AM of 
gas is misaligned with the dark matter then this could poten- 
tially lower the signal. Another example is related to the use 



of weak lensing studies to measure the projected mass den- 
sity of a foreground galaxy in front of background galaxies. 
Since signal from an individual galaxy is weak, to produce 
detectable signals, results of different galaxies are stacked to- 
gether by orienting the images with respect to the shape of the 
central galaxy. If the AM of the galaxies is misaligned with 
respect to the shape of the dark matter halos, then this can 
wash out any ellipticit y signal in the projected mass distribu- 
tions (Bet tetal.l[2QTQh . 

The AM vectors of gas and dark matter, in addition to be- 
ing misaligned with each other, are also not perfectly a ligned 
with them selves within the halo (iBailin & Steinmetj 120051: 
iBett et al.l2010,) . which we refer to as self-misalignment. The 
self-misalignment is fou nd to be most pronounced between 
the inner and outer parts ("B ailin & Steinmetzll2005) . For the 
gas such a self-misalignment could be responsible for warps 
as seen in galactic discs. In recent cos mological hydrodynam- 
ical simulations, Roskar et al.l (120101) show that the warps in 
their disks are due to the misalignment of the AM of the in- 
ner cold gas with that of the outer hot gaseous halo. Hence, 
it is important to understand as to when such a misalignment 
occurs. 

The self-misalignment of AM could also be respon- 
sible for the counter rotating gas as seen in some of 
the galaxies ( Ciri et al.l 119951; i Sil'chenko & MoiseevI 120061: 
ISil'chenkoet a l. 2009). Although recent mergers of gas rich 
systems are generally used to explain them, the models have 
some shortcomings. For example, if the merger is too massive 
it can heat up and thicken the disc considerably; if it is small 
then in some cases it cannot account for all of the counter 
rotating gas (iCiri et al.llT99l iThakar & Rvdenlll996h . Mis- 
aligned AM in galactic halos could provide an explanation 
for this. 

To answer some of the questions posed earlier, we perform 
non-radiative hydrodynamical simulations of merging spheri- 
cal halos and analyze the AMDs of the resulting remnant ha- 
los. We do simulations with various different orbital param- 
eters and study the dependence of the shape parameter a of 
AMDs on these orbital parameters. We also analyze the ra- 
tio Acas/ADM and the misalignment angle 9 of the remnant 
halos. 

An outline of paper is as follows: in Section 2 we describe 
details of setting up initial conditions and methods of extract- 
ing AMDs from halos; in Section 3 we investigate the AM 
properties of these halos. Finally, in Section 6 we summarize 
and discuss our results. 

2. METHODS 

2.1. Initial Conditions and Simulations 

We study binary mergers of spherical halos consisting of 
dark matter and gas. The halos are set up with an exponen- 
tially truncated NFW density profile. 
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Imposing the condition that the logarithmic slope of p at tran- 
sition radius, r = rtr, should be continuous, gives 
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cation gives rise to an extra mass. To compensate for this 
we choose rtr to be less than rvir, such that the total mass of 
the system is mvir- For generating equilibrium realizations 
of the system, comprising of coUisionless particles, we fol- 
low the procedure given by Kazantzidis et al. (2004). In this 
procedure, first the phase space distribution function corre- 
sponding to a given density profile is numerically evaluated 
and then the velocities of the coUisionless particles are as- 
signed by randomly sampling this distribution. The gas is 
setup in hydrostatic equilibrium within the dark matter halo 
assuming a density profile identical to that of the dark matter 

(PGas(^) = P'DM{r)fh/(l - fb), fb = ^baryon /^matter being 

the cosmological baryon fraction). The thermal energy of the 
gas is given by 



u{r) = 



PGas(r) 



^ ,GM{<r) 



(3) 



For all our set ups, we use rd = O.lrvir- Normally, rtr should 
be chosen to be equal to rrvir- However, the exponential trun- 



M{< r) being the cumulative mass enclosed within radius r. 

In Table [T] we list the parameters that are used to set up 1 1 
simulations, each with A^ = 2 x 10^ dark matter particles and 
an equal number of gas particles. Merger parameters were 
selected as follows. A two body merger can be described in 
terms of the motion of a test particle with a reduced mass. A 
bound orbit of such a test particle can be fully characterized 
by semi major axis a and eccentricity e or equivalently by 
orbital time period Torb and spin parameter Xorb of the sys- 
tem. The choice of these parameters is constrained by the fact 
that Torb should be less than the age of universe and that the 
maximum separation between objects rrei = a(l + e) should 
be greater than rviri + ^vir2. The later condition provides a 
lower bound on Torb- For our simulations, given a A we set 
Torb such that rrei = ^viri + ^vir2. The exception being Sim- 
6 which has rrei > ri2 and hence was translated analytically 
till the separation between the halos was equal to ri2. Fur- 
ther details on the setup of merger parameters and its physical 
interpretation is given in Section |A] of the appendix. 

For simulations 1 to 8 we assume the density distribution of 
gas to be same as that of dark matter but in simulations 10 and 
1 1 the gas is allowed to have a different density distribution, 
namely the concentration parameter for gas is different from 
that of dark matter and this is shown in brackets. All the sim- 
ulations except Sim-9 start with non-rotating halos, i.e., zero 
intrinsic spin. For the Sim-9 we use the remnant halo obtained 
from Sim-1 as initial halo and Torb is set to be perpendicular 
to the spin Tint of the halos. The intrinsic halo spins are as- 
sumed to be parallel to each other and are pointing towards 
the z axis. For this setup the direction of orbital AM in spher- 
ical coordinates is given by (^, 6>)orb = (—90,90)°. Three 
other setups similar to this but with (^, O)orh = (—90, 45)°, 
(-90, 135)° and (-90, 180)° were also performed but are not 
listed in Table [TJ 

All simulations were evolved for 10 h~^ Gyr. The simu- 
lations were don e using the sraooth p article hydrodynamics 
code GADGET (ISpringelet al.ll200ll) . By construction, no 
assumptions on a particular background cosmology are made; 
however, for the NFW halo parameters we adopt the concor- 
dance ACDM cosmology with Qx = 0.732, Q^ = 0.267. A 
gravitational softening of 2 kpc h~^ was used. 

In order to compare the AM properties of merger simula- 
tions with those of simulations done in cosmological context 
we additionally use a set of 42 halos (virial masses between 
1.3 X IO^^Mq to 1.5 X lO^^M©), which were selected from 
a 32.5 h~^Mpc box length dark matter simulation (128^ par- 



TABLE 1 
The initial set up parameters of merger simulations and their final properties. 



Sim 


Mtot 


fm 


Aorb 


^initial 


^orb 


Aint 


mvir 


/;; 


Cfinal 





Agas 


Adm 


Agas/AoM 


«DM 


«Gas 


«Gas/«DM 


1 


100 


0.5 


0.05 


10.0 


6.70 


0.0 


83.6 


1.0 


10.4 


0.4 


0.037 


0.039 


0.96(1.17) 


0.85 


0.97 


1.14 


2 


100 


0.5 


0.05 


5.0 


6.70 


0.0 


83.4 


1.02 


5.75 


0.6 


0.042 


0.039 


1.08(1.27) 


0.87 


1.08 


1.24 


3 


100 


0.5 


0.05 


15.0 


6.68 


0.0 


84.4 


1.0 


15.1 


0.8 


0.035 


0.038 


0.91(1.10) 


0.78 


0.90 


1.15 


4 


100 


0.5 


0.01 


10.0 


6.40 


0.0 


84.2 


1.0 


10.8 


1.7 


0.0079 


0.0084 


0.94(1.04) 


0.89 


0.93 


1.04 


5 


100 


0.5 


0.10 


10.0 


7.50 


0.0 


83.3 


1.03 


11.3 


0.4 


0.093 


0.074 


1.25(1.30) 


0.86 


1.02 


1.19 


6 


100 


0.5 


0.05 


10.0 


10.0 


0.0 


80.2 


1.02 


10.5 


0.1 


0.043 


0.037 


1.15(1.49) 


0.79 


0.96 


1.22 


7 


100 


0.1 


0.05 


10.0 


7.40 


0.0 


91.8 


0.95 


9.2 


2.0 


0.046 


0.024 


1.96(1.62) 


0.82 


0.74 


0.90 


8 


100 


0.3 


0.05 


10.0 


6.60 


0.0 


86.6 


0.99 


10.5 


1.2 


0.044 


0.037 


1.17(1.28) 


0.85 


0.86 


1.01 


9 


167.2 


0.5 


0.05 


10.08 


6.70 


0.039 


137.5 


1.04 


11.5 


18.2 


0.044 


0.041 


1.09(1.22) 


0.75 


0.94 


1.25 


10 


100 


0.5 


0.05 


10(1) 


6.70 


0.0 


79.3 


0.80 


8.9 


2.2 


0.029 


0.036 


0.79(0.68) 


0.82 


1.02 


1.24 


11 


100 


0.5 


0.05 


10(25) 


6.70 


0.0 


84.5 


1.06 


14.3 


1.5 


0.036 


0.041 


0.88(0.94) 


0.81 


0.91 


1.12 



Note. — In the table, columns 2 to 7 are the parameters used to set up the initial conditions of the merger simulations. The columns are as follows: the total 



mass of the simulated system Mtot = rni + m2 (10-*^^ h ^ Mq), the fractional mass of the least massive halo fn 



niin(mi ,rn2) 
(mi+m2) 



, the initial spin parameter 



of the whole system Aorb» the concentration parameter of the initial halos Cinitiah the orbital time period of the merging system Torb ( h -"^ Gyr), and the 
intrinsic spin of merging halos Aint • Note, for Sim- 10 and Sim- 11 the Cinitiai for gas is different from that of dark matter, hence, the Cinitiai for g^s is quoted 
in parenthesis. Columns 8 to 17 describe the properties of the merger remnants formed by the simulations. Here mvir (10"'^'-' h"-*^ M©) is the virial mass of the 
remnant , /^ is the baryon fraction within virial radius relative to cosmological baryon fraction, Cfinai is the concentration parameter of the remnant, 6 (degree) 
the misalignment angle between the dark matter and gas AM vectors. These are followed by the spin parameter A and the shape parameter a of AMDs for gas 
and dark matter. The quantity in brackets is the spin ratio at t = 6 h~^ Gyr. The spin parameter is calculated using the definition given bv ' Bullock et a l. (2001). 



tides), and were resimulated with gas at higher resolution by 
ISharm a & Ste inmetzl (l2Q05h using GADGET. In these halos 
the number of dark matter particles within the virial radius 
ranges from 8000 to 80,000. A gravitational softening of 
2 kpc h~^ was used. 

2.2. Calculation of Angular Momentum Distributions 

Dark matter particles are assumed to be collisionless and 
thus a significant amount of random motions are superim- 
posed onto the underlying rotational motion. So in order 
to calculate AMDs , the v elocity has to be smoothed (see 
ISharma & Steinmet3l2005|) . Since the rotational motion is 
very small compared to the random motion, one needs to 
smooth with a large number of neighbors. This large scale 
smoothing introduces systematic biases which needs to be 
taken into account. Smoothing the Cartesian components of 
velocity spuriously underestimates the rotation for particles 
near the axis, as < Vx >=< V y >=< Vy >^ ne ar the 
axis. To avoid this problem in Sharma & Steinmetz (2005) we 
smoothed the Cartesian components of AM vectors instead of 
velocities. As we will demonstrate later, the angular velocity 
^ is nearly constant near the center. This implies that the AM 
vector j has a strong, monotonically increasing radial depen- 
dence on cylindrical coordinate Tc. This results in an over- 
estimate of the AM of particles close to the axis. Existence 
of a strong radial density gradient further leads to underesti- 
mate of AM for particles along the equator. To reduce some 
of these problems, in this paper we choose to smooth the z 
component of the angular velocity vector 0,^ (the halo being 
oriented such that the z axis points along the total AM of the 
component being smoothed). A simple top hat kernel is used 
for smoothing. The number of smoothing neighbors was cho- 
sen to be 5 X 10~^ times the total number of particles within 
the virial region. This makes the smoothing volume indepen- 
dent of the number of particles in the halo. Note, smoothing 
is only employed to calculate the shape parameter a of the 
resulting AMDs. 

3. ANGULAR MOMENTUM PROPERTIES OF REMNANT HALOS 

The final properties of the merger remnants are given in Ta- 
ble [T] We note that the virial mass mvir of the remnant is less 
than the total mass of the system Mtot • Hence, a fraction of 



mass is lost. The lost mass fraction is an increasing function 
of the kinetic energy KE involved in the collision and a de- 
creasing function of the total binding energy of the system. A 
detailed description of the mass structure of remnant halos is 
given in Section |B] in appendix. In the present section we ad- 
dress the AM properties of the remnant halos. First, we study 
the angular velocity and AMDs of the merger remnants and 
compare them with those of halos simulated in cosmological 
context. Next, we study the evolution of AM with time of our 
fiducial equal mass merger, Sim-1, and shed light on the phys- 
ical mechanisms responsible for distributing the AM within 
the halo. This is followed by studying the dependence of AM 
properties on orbital parameters. In the penultimate subsec- 
tion we look into the issue of misalignment of AM vectors 
and in the last subsection we study the spatial distribution of 
low AM material. 

3.1. Angular Velocity and Angular Momentum Distribution 

of Halos 

We explore the angular velocity 1] as a function of spherical 
coordinates r and 9 for the remnant halos at t = 10 h~^ Gyr, 
which represents the final relaxed configuration. By angular 
velocity we mean the z component of AM, with z axis point- 
ing along the direction of the total AM vector of the compo- 
nent being analyzed (gas or dark matter). The angular velocity 
Vt of both gas and dark matter is found to be nearly indepen- 
dent of 9, both for r < rvir and r < rvir/2 (lower two panels 
of Figure [TJ. The weak trend that exists is monotonic and 
favors faster rotation towards the equator. This suggests that 
shells of matter are in solid body rotation. The top panel in 
Figure [T] shows the radial profiles of gas and dark matter. In 
general 1] is a decreasing function of radius r but the profiles 
seem to flatten for r < 0.2rvir. As compared to dark mat- 
ter, the gas is found to rotate faster in the inner regions and 
slower in the outer regions. For comparison the angular ve- 
locity profiles of halos simulated in cosmological simulations 
(from Sharma & Steinmetz (2005)) are shown in Figure |2] As 
in merger simulations they are nearly independent of angle 
6>, are a decreasing function of radius r, and show faster ro- 
tation for gas in the inner regions. However, the faster ro- 
tation for gas is not as strong as in merger simulations and 
the dip in gas rotation at about r = O.Trvir is also not seen. 
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Fig. 1. — Median angular velocity Q = jz/ix'^ +2/^) as function of 
radius r and angle 6 for halos formed by mergers (merger simulations 1 to 
8, excluding 4). ^2 as a function of r and was calculated by binning the 
particles so as to 1000 particles per bin. The dashed lines show 16th and 
84th percentile values. The angular velocity profiles seems to flatten out for 
r < 0.2r/rvir • In the top panel it can be seen that in the inner regions the gas 
rotates faster than dark matter. Note, the gas profiles are much more smooth 
than that of dark matter and this is because the dark matter has significant 
amount of random motion superimposed on the actual rotation which is quite 
small. 



Note, these are median profiles: on a one to one basis the gas 
and dark matter can show much more prominent differences 
as is revealed by the fact that there is significant scatter in 
the ratio ^Cas/^DM- Also a real halo has much more com- 
plex merger history which can probably reduce the difference 
between dark matter and gas in the inner regions. The dif- 
ference in the outer parts is probably due to the fact that the 
initial merging halos have an exponential cut off in the outer 
parts whereas in real simulations the halos are much more ex- 
tended and moreover there is also smooth accretion onto the 
halos. Hence, the outer parts of merger remnants may not be 
an accurate representation of the real halos. 

We now study the AMDs of the remnant halos at t = 
10 h~^ Gyr. The AM of each particle is obtained by smooth- 
ing its angular velocity with 400 neighbors. For fitting the 
A MDs we use the following a nalytical function (for details 
see lSharma & Steinmet"3l2005h 



P{J)- 



1 



Jd^^) 



(j)" le ^/^'^where jrf= jtot/c^. (4) 



jtot being the mean specific AM of the system. Writing P in 
terms of s = j/jtot and replacing jd the cumulative distribu- 
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Fig. 2. — Median angular velocity Q = jz/(x'^ +2/^) as function of 
radius r and angle ^ for 21 halos simulated in a cosmological context. The 
halos where selected so as to have more than 30,000 particles within the virial 
radius individually for both gas and dark matter. ^2 as a function of r and 6 
was calculated by binning the particles so as to have 1000 particles per bin. 
The dashed lines show 16th and 84th percentile values. The angular velocity 
profiles seems to flatten out for r < 0.2r/rvir- In the top panel it can be 
seen that the gas rotates faster in the inner regions than that of dark matter. 



tion reads as 



P(< s) =7(0^, as) 



(5) 



where 7 is the I ncomp lete Gamma function. In 

ISharma & Steinmet3 (|2005h this function was used to fit the 
AMDs of halos obtained in cosmological simulations and of 
model exponential disks embedded in NFW potentials. It was 
found that for AMDs of exponential disks embedded in NFW 
potentials the shape parameter a is greater than 1.3 whereas 
for cosmological halos values are typically smaller than 1 
(< Q^DM >= 0.83 and < Ofgas >= 0.89 ). For dark mat- 
ter for fiducial Sim-1 we find a = 0.85 whereas for others 
it is given by 0.75 < a < 0.9. For gas in Sim-1, a is 0.97 
and for others it is between 0.74 and 1.08. The gas has sig- 
nificantly larger a than dark matter and this is because of the 
fact that the gas rotates faster than dark matter in the inner 
regions. Merger simulations successfully reproduce the fact 
that acas > <^DM as in cosmological simulations. If we take 
Sim-1 as the fiducial case then for dark matter the value of a 
is in excellent agreement with cosmological simulations but 
for gas we find that it is about 8% higher. As discussed earlier 
the gas in merger simulations are an idealized case and in real 
halos the gas rotation profiles are slightly flatter in the outer 
parts and this explains the slightly lower a in them. 
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Fig. 3. — The top curve shows the evolution of the half mass radius 
^median- The bottom three curves show the evolution of the specific angu- 
lar momentum j of the gas and dark matter components for Sim-1 (an equal 
mass merger) for the full system and separated into the inner half mass and 
outer half mass. 



3.2. Evolution of Angular Momentum with Time 

In this section we analyze the time evolution of the specific 
AM of gas and dark matter components for Sim-1 (Figure |3}, 
which is an equal mass merger of halos with concentration 
parameter of 10, Torb = 6.7 h~^ Gyr and Ainitiai = 0.05. We 
use this as a fiducial case to understand the main properties 
of the evolution, specifically the origin of the differences in 
the AM properties of gas and dark matter. We analyze the 
evolution of the total AM as well as that of the inner and the 
outer parts of equal mass, separated by median radial distance. 
The evolution of median radius Tmedian, separating the two 
halves is also plotted alongside. We divide the evolution into 
four stages. Stage 1 from — 2 h~^ Gyr, Stage 2 from 2 — 
3 h~^ Gyr, Stage 3 from 3 — 6 h~^ Gyr and Stage 4 from 
6-lOh-i Gyr. 

During a merger the system first collapses to a compact 
configuration (Stage 1 and Stage 2), marked by a decrease 
in rmedian (scc also Figure |4]). In Stage 3 the system expands, 
as shown by slight increase in rmedian, and then in Stage 4 the 
system evolves without any significant change in the density 
structure. It can be seen from the evolution of the total AM of 
the system (middle lines in Figure [3]) that in Stage 2 the gas 
loses about 10% of its AM to dark matter. The dark matter 
gains AM in this stage but its change is quite small since the 
mass of dark matter is much larger than that of gas. In Stage 3 
and 4 the total AM of gas and dark matter is nearly constant. 
In contrast, when the AM of inner and outer half masses are 



analyzed separately, significant differences can be seen. For 
dark matter there is an inside-out transfer of AM in stages 1 , 2 
and 3. Due to dynamical friction the inner part loses AM con- 
tinuously to the outer part until it virializes to form a pseudo 
equilibrium distribution after which the evolution stops. It can 
be seen that the AM evolution of gas is decoupled from dark 
matter, from stage 2 onwards. Initially, both the inner and 
outer parts of gas lose AM to dark matter. However, in Stage 
3, when the inner parts start to expand, for the gas the AM 
is transferred to the inner parts from the outer parts. The fact 
that the rise in AM of gas in inner parts is almost the same as 
the fall in AM of gas in outer parts means that the transfer of 
AM is purely between the gas components. This transfer is 
because the expanding inner part of gas that also has low AM, 
shocks with in-falling outer part that has high AM, thereby 
leading to transfer of AM. This is visible more clearly in Fig- 
ured] where we plot the velocity field in the X-Y plane within 
a radius of 125 kpc and \z\ < 20. At 3.0 h^^ Gyr the halos 
can be seen crossing each other and at 3.2 h~^ Gyr they have 
crossed and are now pushing against the outer material of the 
other halo which is still falling in. The outer material falling 
in from upper right and lower left comers pushes and trans- 
fers AM to the expanding inner regions. With time the shocks 
progressively move outwards. 

In contrast the dark matter cannot shock, their particles can 
cross each other and they exchange energy and AM via violent 
relaxation. It is easier to understand their evolution in terms of 
an inside-out spherical collapse simulation in which the inner 
regions collapse faster than the outer regions. In such a sys- 
tem as described in B inney & Tremain e (2008) a high energy 
particle in outer region falls into a gradually steepening po- 
tential and hence gains kinetic energy. Later when it starts to 
move out the inner region has already expanded and hence it 
has to climb out of a shallower potential. The net result of all 
this is that a high energy particle that falls in late gains energy. 
Now, the impact parameter of the particle during collision is 
high for particles in the outskirts that are falling in late. Since 
the AM of a particle is proportional to the impact parameter 
it is also high for them. Hence, high AM particles mostly end 
up orbiting in the outer regions. In contrast for gas these late 
falling high AM material shocks and transfers its AM to the 
inner parts. 

Overall the conclusion is that due to gas dynamical effects 
the baryons are more efficient in depositing the AM to the 
inner parts of the halo making it rotate faster than dark matter 
in the inner regions. The faster rotation of the gas in the inner 
region is visible in Figure [TJ and also in the bottom panels of 
Figure m The faster rotation of gas also responsible for acas 
being greater than q^dm- 

The above hypothesis suggests that increasing the energy 
of the collision or the orbital AM should make the outside- 
in transfer of AM for gas and inside-out transfer of AM for 
dark matter more stronger. The Sim-6, which is same as Sim- 
1 except for the fact that it has higher Torb meaning more 
energetic merger, does reveal this. As expected, the ratios of 
AGas/^DM and aGas/<^DM are found to be higher for Sim-6. 
Similarly, when considering Sim-4, Sim-1 and Sim-5 having 
lowest intermediate and highest orbital AM one finds that the 
ratios of Acas/ADM and acas/o^DM increases monotonically 
(see Table [TJ. These results thus provide support to the above 
hypothesis. 

In the final Stage 4 of the evolution the system has almost 
reached a pseudo equilibrium. During this stage, for the gas 
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Fig. 4. — Velocity vector field in the X-Y plane along with density map for gas and dark matter at various stages during the evolution of an equal mass merger 
simulation (Sim-1). Field is shown for particles within a radius of 125 kpc and \z\ < 20. The maximum length of the arrow corresponds to 250 km s~^ and is 
shown in left hand comer of each panel. At 3.2 and 3.4 h~^ Gyr it can be seen that at the regions near the shocks (top right and lower left) the in-falling gas in 
the outer parts is transferring angular momentum to the inner parts. By t = 6 h.~^ Gyr the gas can be seen to rotate faster in the inner regions as compared to 
dark matter. 

there is a gradual transfer of AM from the fast rotating inner 
layers to the slow rotating outer layers. 

3.3. Dependence of Spin Ratio Acas/ADM on Orbital 
Parameters and its Evolution with Time 

Having understood the AM evolution in the inner and outer 
parts, we now try to understand the evolution of AM within 
the virial radius, which is commonly employed to measure the 
spin of the halos. Figure [5] describes the evolution of the spin 
ratio Acas/ADM with time for various different merging sce- 
narios. Note, the spin parameter is calculated using the defini- 
tion A = jtot/V^Gmvir^vir (.Bullock ct al.ll2QQlh . where jtot 



is the specific AM of the material within virial radius rvir, 
and mvir is the virial mass. At each stage of the evolution 
we identify the virial region by means of the spherical over- 
density criterion and then compute the relevant properties of 
the virialized remnant halo. In Stage 1 (0 — 2 h~^ Gyr) of the 
evolution the ratio is close to 1. In Stage 2 (2 — 3 h~^ Gyr) the 
ratio drops by about 10 - 20%. In Stage 3 (3 - Gh"^ Gyr), the 
ratio rises and reaches a peak at around 6 — 7h~ ^ Gyr and then 
in Stage 4 the ratio decreases (except for Sim-7). It is easy to 
understand the time evolution of spin ratios in the context of 
the discussion done earlier in Section 13.21 In Stage 1 the gas 
and dark matter have not yet decoupled so the ratio is close to 
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Fig. 5. — Evolution of Acas/^DM with time for various merging scenar- 
ios. 

1. In Stage 2 the gas loses its AM to dark matter and hence a 
drop in the spin ratio. In Stage 3, in the inner regions the dark 
matter loses AM while the gas gains, this results in a rise in 
the spin ratio. Finally, in Stage 4 the AM of dark matter in the 
inner regions remains nearly constant whereas for gas there 
is an inside-out transfer and this again results in a drop in the 
spin ratio. 

Next, we study the dependence of spin ratio on the orbital 
parameters. In each of the panels in Figured] we vary one of 
the orbital parameters (namely /m, c, A and Torb) while keep- 
ing the other parameters identical to that of benchmark Sim- 1 . 
In top left panel we compare Sim-1, Sim- 8 and Sim-7, having 
fm =0.5, 0.3 and 0.1 respectively, fm = 1712/ {mi + 777.2) be- 
ing the mass fraction of the smaller merging halo. At a given 
time the gas to dark matter spin ratio is found to be higher for 
a lower value of fm- For fm = 0.1 it continues to increase 
even in Stage 4 and reaches a value as high as 2. 

In the second panel, i.e., top right, we plot the results 
for mergers with different values of concentration parameter, 
Sim-1, Sim-2 and Sim-3 having cinitiai =10.0, 5.0 and 15.0 
respectively. Lower concentrations yield higher spin ratios. 
It can be seen from Table [T] that Adm is largely unaffected by 
the change in cinitiai whereas Acas increases with lowering the 
concentration. In Sim- 10 and Sim- 11 we vary the concentra- 
tion parameter of gas, setting it to 1 and 25 respectively, and 
keep the concentration of dark matter constant at 10.0. The 
Sim- 10 is designed to mimic the case of a halo where the gas 
is puffed up by feedback from star formation whereas Sim-1 1 
mimics the case where the gas has cooled and collapsed to the 
central regions. Table [T] shows that when considering the to- 
tal AM content, the concentrated gas loses more AM than the 
puffed gas. This demonstrates the AM catastrophe problem in 
which due to excessive cooling the gas gets concentrated and 
during subsequent evolution lose AM as a result of dynamical 
friction. Surprisingly, when AM is measured with in the virial 
region the puffed up gas has less AM. This is because for the 
puffed up case significant amount of gas is outside the virial 
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Fig. 6. — Evolution of aoas with time for various merging scenarios, a is 
calculated from the AMDs obtained from smoothing the motion of particles 
obtained from simulations. 

radius and this gas also has high AM whereas for the con- 
centrated gas case all the gas ends up within the virial radius. 
This is reflected in the baryon fraction as shown in Table [T] 
which is 0.8 for the former and 1.06 for the later. 

In the bottom left panel we look at the role of varying the 
orbital AM. We compare Sim-1, Sim-4 and Sim-5 having a 
Aorb =0.05, 0.01 and 0.10 respectively. Increasing Aorb be- 
yond 0.05 increases the spin ratio while lowering it does not 
affect the results significantly. Finally, we investigate the role 
played by the kinetic energy associated with the collision, 
which is controlled by varying the parameter Torb- Larger 
Torb means that halos approach each other from a farther dis- 
tance and have more energetic collision. We compare Sim-1 
and Sim-6 which have Torb = 6.7 and 10.0 respectively. More 
energetic collision leads to higher spin ratio. In light of dis- 
cussion in Section [T2I the effects discussed above are due to 
the fact that late infalling gas, in case of high Torb and Xorb, 
shocks more strongly leading to more transfer of AM to inner 
parts and for dark matter the late infalling gas is more ener- 
getic and is more likely to escape outside the virial radius. 

In general it can be seen that at around 6 h~^ Gyr i.e., 
3 h~^ Gyr after the merger the ratio AGas/^DM > 1 for all 
merging scenarios and this provides an e xplanation for the 
results of ISharma & Steinmetzl (l2005h and lChen et"an (l2003h 
where they find AGas/^DM ~ 1.4 and 1.2 respectively for 
halos simulated in a cosmological context. 

3.4. Dependence of Shape Parameter a on Orbital 
Parameters 

In this section we explore the role of the orbital parameters 
on the shape parameter a of gas obtained by fitting the AMD 
of the remnant halos by Equation (O. In Figure |6] we plot 
the evolution of the shape parameter a Gas for various merg- 
ing scenarios. The comparisons done in various panels are the 
same as in FigureO The values of o^Gas below t = 3 h~^ Gyr 
are not relevant for the study here since the merger has not yet 



happened. Between 3 — 5 h~^ Gyr there is a sHght variation 
where the halo is still relaxing, but beyond that for all orbital 
geometries a has very little evolution with time (Figure [6]). 
As an apparent trend, a decreases slightly with time ( except 
Sim-2). Varying the parameter Aorb or torb does not seem to 
affect the values of a. Decreasing the mass ratio fm decreases 
the value of a, while decreasing the concentration parameter 
Cinitiai increases its value. Varying only the concentration of 
gas as in Sim- 10 and Sim- 11 also has similar effect (see Ta- 
ble[T]), namely puffed up halos have higher a whereas concen- 
trated halos have lower a. In the context of the AMD problem 
this means puffing up gas by means of feedback can partially 
help to resolve the problem, but the value of a = 1.02 is still 
far short of that required to form exponential disks (a > 1.3). 
Hence, just by itself the puffing up of gas is not enough to 
solve the problem. 

Finally, we note that the gas in general has higher a than 
that of dark matter and the aGas/<^DM increases with in- 
crea se o f Torb and Aorb (see Table [TJ. As discussed in Sec- 
tion 13.21 this is due to the differences between gas dynamics 
and collisionless dynamics. The fact t hat ac.?^s is greater than 
QfDM, is consistent with the findings of 'Chen et al.l (|2003l) and 
[Sharma & Steinmetz (2005) for cosmological halos, and our 
results here provide an explanation for it. 

3.5. Misalignment of the angular momentum vectors of gas 
and dark matter 

The misalignment angle 6 for all simulations is tabulated in 
Table [T] Mergers with zero intrinsic spins do not seem to gen- 
erate any significant misalignment in the final remnant halo. 
The misalignment angle is less than 2° for all orbital geome- 
tries. Now, if we imagine a merger of halos with non-zero in- 
trinsic spins, then the final angular momentum can be though 
of as j = jint + jorb, jint being the intrinsic specific AM and 
jorb the orbital specific AM. Due to differences between gas 
and collisionless dynamics we already know j^^® > j^^, 
similarly the contribution of the intrinsic components could 
also differ. If jint is aligned with jorb then there is no possi- 
bility to generate any misalignment. However, if they are not 
aligned then we can expect to see misalignment, except for 

the special case where jSbVj?r^ = J?nT/J?X- 

To test the above scenario, in Sim-9 we merge two halos 
(extracted from Sim-1) having non-zero spin and an orbital 
AM which is perpendicular to the spin. The remnant halo 
is found to be significantly misaligned with a misalignment 
angle close to 17°. In Figure [7] we show the orientation of 
the AM vectors of gas and dark matter as measured in radial 
shells for various different directions of the orbital AM. We 
mainly concentrate on regions with r > O.lrvir which should 
be quite reliable given that our gravitation softening is about 
O.Olrvir- The intrinsic spin has the direction 6 = 0.0. The 
solid and dashed lines are the differential profiles while the 
rest are cumulative profiles. In the top panel the lower hori- 
zontal line marks the mean expected 6 for the halo assuming 
uniform mixing. The upper line shows the angle for the or- 
bital AM. The gas and dark matter show very different trends. 
For dark matter the inner region is dominated by orbital AM, 
the outer by intrinsic spin and the middle region has interme- 
diate direction. For the gas the inner region has intermediate 
values, the outer region is dominated by orbital AM and the 
middle region is dominated by intrinsic spin, which points to- 
wards 6> = 0°. In the rightmost column corresponding to a 
retrograde merger the gas even shows a spin flip in the middle 



regions. 

In the bottom panels it can be seen that the cumulative mis- 
alignment angle defined as /3DM-Gas(< r) = cos~^(Jdm(< 
^)Jdm(< r)) increases inwards into the halo. The cumula- 
tive self-misalignment of the AM , /3dm-dm and /^Gas-Gas, 
which is measured with respect to the AM within the virial 
radius also shows similar trend. In the bottom row the mag- 
nitude of misalignment increases from left to right, i.e., with 
increase of angle between orbital and intrinsic AM. 

In the panels in second column, the case of 6>orb = 90, it can 
be seen that most of the misalignment is due to the (j) direction 
varying sharply in the inner regions. Moreover, the gas and 
and dark matter AM vectors seem to be pointing in opposite 
directions in (j). This is surprising given the expected value of 
(j) is —90°. It is not clear if the gas is being torqued by dark 
matter or is it simply rearrangement of AM. To understand 
this cumulative profiles are plotted as dotted lines. By r = 
Tvir the seems to have averaged to the expected value for 
both gas and dark matter. Also, in 6 the cumulative profiles 
tend to the expected value at large r for both dark matter and 
gas. This suggests redistribution and self torquing to be the 
main mechanism for the variation of the gas AM direction. 
However, beyond rvir, ^ for gas is slightly larger than 90, 
hence some amount of torque must have been exerted on it 
from dark matter. The panels in other columns also lead to 
similar conclusion. 

We now compare the results of cosmological simulations 
with that of merger simulations. In Figure [8] in the top panel 
we plot the cumulative misalignment angles /3dm-dm(< 
^), /^Gas-Gas (< r) and /3DM-Gas(< r) as defined earlier. The 
trends in the top panel are similar to the trends in the second 
column of bottom row in Figure |7] which corresponds to the 
most probable orientation of a merger. Figure |8] shows that 
the median misalignment of gas with respect to dark matter is 
about 20° which is reproduced by Sim-9. Note, higher mis- 
alignments can also be achieved if Oorh is greater than 90°. 
The fact that results of the cosmological simulations are suc- 
cessfully reproduced by the merger simulations leads us to 
conclude that the difference in gas and collisionless dynam- 
ics is the main cause of misalignment of AM vectors as seen 
in cosmological simulations. Moreover, misalignments occur 
only when the intrinsic spins are not aligned with the orbital 
AM. 

Top panel in Figure [8] also shows that the misalignments are 
more pronounced in the inner parts than the outer parts. The 
fact that the total AM vectors are dominated by the AM in the 
outer parts is partly responsible for this. Finally, in the ori- 
entation profiles of cosmological simulations the dark matter 
shows more misalignment than the gas whereas the opposite 
was true for the merger simulations. Given that a real halo 
has a much more complex merger history than that of a single 
merger as shown here, we do not consider the discrepancy to 
be too significant. Moreover there is a significant scatter about 
the median profiles as shown in Figure [8] which means that on 
a one to one basis the gas and dark matter can have different 
trends suggesting that they are sensitive to the merger history 
and hence could be employed to understand them. 

We now compare the amount of misalignment seen in our 
cosmological halos with those of previous studies. The pre- 
vious studies were mostly confined to dark matter only sim- 
ulation s hence only /3dm -dm can be compared. iBett et al.] 
(1201 Oh show the cumulative misalignment angle with re- 
spect to AM vector of material within r < 0.25rvir while 
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Fig. 7. — Orientation angles 6 and cj) of angular momentum vectors (of gas and dark matter) for different merging scenarios. In each case the intrinsic spin of 
merging halos points along the z axis. Four cases with different directions of orbital angular momentum are shown. In the upper two panels, the horizontal black 
solid line shows the direction of the total angular momentum of the system while the dot dashed line shows the direction of orbital angular momentum. The curves 
were obtained by binning the particles in 25 logarithmically spaced bins in radial distance. The dotted lines show the cumulative profiles while the solid lines are 
differential profiles. The dashed line is the differential profile of the whole system (gas and dark matter). The bottom row shows the cumulative misalignment 
angle of dark matter with itself, gas with itself and between gas and dark matter. The gas and dark matter angular momentum vectors are significantly misaligned 
with each other, especially in the inner regions. Additionally, in the inner region r < 0.2rvir, the angular momentum vectors of gas and dark matter are also 
individually misaligned with their total angular momentum vectors. 

formation and feedback could further alter the angular mo- 
mentum of the disc. Hence, in genera l we expect real galaxies 
to be even more misaligned. In fact, iBett et al.] (12010 ) using 
simulations with star formation and feedback , find the median 
misalignment between the AM of the central galaxies and the 
DM halo to be -30°. 



iBailin & Steinmetj (120051) use a differential distribution. In 
order to facilitate comparison with earlier studies we addi- 
tionally plot in middle and bottom panels the cumulative 
/3{> r) profiles and the differential /3(r) profiles. For the 
range of masses considered here Bett et al. (2010) find a value 
of around /3dm-dm(^ < 0.25rvir) = 25° (their Fig 4). 
The top panel of our figure shows the corresponding quan- 
tity to be 30° which is in good agreement with their results. 
iBailin & Steinmetj (1200 5) plot the orientation profiles with 
respect to AM measured in different shells. If the shell at 
r = 0.6rvir is taken be the representative of the total AM, 
then this gives a value of 25 and 35 for /3DM-DM(0.1rvir) 
and /3DM-DM(^vir) which is again in very good agreement 
with profile shown in bottom panel of our Figure [8] 

Our result for misalignment between gas and dark mat- 
ter, median /3DM-Gas(<liw) =^°, are in agreement with 
those of van den Bosch et al.l (120021) . This was also reported 
in ISharma & Steinmeta (|2005|) ! These results are for non- 
radiative gas. A portion of this gas would cool and later on 
form disc galaxies. Additional, physical processes like star 



3.6. Spatial Distribution of Low and Negative Angular 
Momentum Material 

The AMD of gas in halos simulated in a cosmological con- 
text, shows an excess of low AM material as compared to the 
AMD required to form an exponential disc. If 5 = j/jtot 
is the specific AM normalized to the mean specific AM, then 
< s < 0.1 is the typical region where the theoretical predic- 
tion differs from that of the exponential disc. Hence, we select 
particles in this range and study their distribution in space. In 
addition, there is also the issue of material with negative AM, 
which can arise from two sources. The first source is random 
turbulent motions and the second source is large scale flows 
that are remnants of shocks and misaligned minor mergers oc- 
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Fig. 8. — Angular momentum orientation profiles of gas and dark matter 
for halos simulated in cosmological context. The top panel shows the cu- 
mulative profile with /3(< r) = j(< r).jvir. the middle also shows the 
cumulative profile but for /3(> r) = j(> r)-jvir and the lower panel shows 
the differential profile with /3(r) = j(^)-jvir- The angle between angular 
momentum vectors of gas and dark matter is defined as /3dm -gas (< ^) = 
Jdm(< r').JGas(< r). ^DM-gas(> r) also having an analogous definition. 



curring after the major merger. The negative AM due to the 
former source would be typically in regions with low AM and 
would vanish when velocities are smoothed locally, as is done 
while calculating the AMDs. On the other hand, large scale 
flows cannot be ea sily smoothed and is the main reason why 
ISharma & Steinmetz (2005) find that, in spite of smoothing, 
cosmological halos have about 8% of matter in negative AM. 
During the assembly of the disk the negative AM is going to 
further enhance the fraction of low AM material, hence it is 
also important to study its distribution in the current context. 
In Figure |9] we show the x-z and y-z density maps of low 
and negative AM gas particles as defined above in various ha- 
los. AM is computed from raw un-smoothed velocities. In 
the plots the z axis is aligned with the total AM vector of 
the halo. The top two rows are for halos from cosmological 
simulations whereas the lower two rows are for merger sim- 
ulations. Among these, the third row is for a merger where 
intrinsic spins are misaligned with the orbital AM (Sim-9), 
and the fourth row is for the fiducial case of an equal mass 
zero intrinsic spin merger (Sim-1). The plots show that the 
low AM material is near the center and in a conical region 
around the rotation axis. The negative AM material is also 
in regions where the AM is low. As discussed earlier in Sec- 
tion 13.11 angular velocity ft is nearly independent of angle 6 



and is a decreasing function of radial distance r. In terms of 
cylindrical radius Re = \/ x^ + y^ the AM is given by l^i?^. 
Hence, the AM is low along the axis of rotation. The conical 
shape is due to the fact that ^^ is a decreasing function of r. 
Assuming negative AM is due to random turbulent motions, 
one expects it to be in regions where ^i^/^random is small, 
which would again be similar to the distribution of low AM 
material. 

For the remnant halo in Sim-1, a part of the negative AM 
material of gas is distributed in a ring shaped structure in 
the x-z plane. During a merger, a plane of compressed and 
shocked gas is formed which is ejected out radially. The ring 
is created when such a gas which has very low AM falls back 
at a later time. Note, only about 3% of the gas is in such form 
which after smoothing reduces to 1%. 

The spatial distribution of low AM material has a depen- 
dence on the merger history of the halos. For example, in the 
third row, which is a merger of halos with intrinsic spins mis- 
aligned with the orbital AM, the central region looks more 
puffed up in x-z projection in comparison to the halo in the 
fourth row. In y-z projection one can see the that the central 
region is twisted. This appearance is because the AM in the 
inner regions is misaligned with respect to the total AM. The 
cosmological halo in the top row also shows such a behavior 
suggesting a major merger with misaligned spins. The halo 
in the second row also has slightly twisted axes in the inner 
region but is very similar to the halo in the bottom row sug- 
gesting that the intrinsic spins of its progenitors were either 
small or well aligned with the orbital AM. 

The characteristic distribution of low AM material found in 
remnant halos as well as cosmological halos suggests that dur- 
ing galaxy formation a mechanism which preferential ejects 
material from the central regions and prevents further mate- 
rial from collapsing along the rotation axis may alleviate the 
AMD problem. Such preferential ejection may be possible 
with feedback from star formation or AGN. Essentially, the 
inner parts would collapse first and start forming stars. The 
feedback would then drive a radial outflow, but since the as- 
sembling gas will have a flattened configuration with density 
being highest near the equatorial regions, the outflow would 
naturally be stronger along the poles causing preferential ejec- 
tion of low AM material. 

To further investigate this idea, we devise a simple geomet- 
ric criteria to selectively remove the low AM material and then 
check the resulting AMDs to see if they conform to those of 
exponential discs, i.e., a > 1.3. Firstly, we fit the the angular 
velocity as a function of radial distance for both gas and dark 
matter by a simple analytical expression. The profiles were 
found to be well approximated by the following equation (see 
Figure fTOl). 
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Here, Avir = IStt^ + 82(1]^ - 1) - 39(l]m - 1)^ is the virial 
over-density parameter used to calculate th e virial region for 
a given matter density Q^ of the universe (iBryan & NormanI 
1199 8h . and H = 0.1 kms"^kpc"^h is the bubble constant. 
Note, jf'tot is the mean specific AM of each component. For 
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Fig. 9. — Spatial density maps of particles with low AM, i.e., < s < 0.1 (columns one and two) and negative AM, i.e.,s < 0) (columns three and four) in 
various halos within the virial radius. Particles are shown in x — z and y — z plane with z axis pointing in the direction of angular momentum. The fraction of 
low and negative AM particles is also labelled on each plot. The grey scale showing the density maps is normalized to the maximum density in each plot. The 
top two panels are for halos from cosmological simulations while the lower two panels are for halos formed by merger simulations, namely, Sim-9(third row) and 
Sim-1 (fourth row). The merger simulations results are shown for the final relaxed configuration at t = lOh"-*^ Gyr. Particles with low and negative angular 
momentum are concentrated in the center and along the axis of rotation. 
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0.267 the expected specific AM is then given by 

s = - — = 35.61- — - 

Jtot 1 + [r/ro) 



1.75 



(8) 



where Re = \/x^~+~^ is the cyHndrical radius. Particles 
with s' > So are selected for removal. The resulting ge- 
ometry for various values of sq is shown in Figure \TT\ The 
geometry is conical in shape and resembles those seen in out- 
flows. For a fiducial setting of 5o = 0.5, the AMDs before 
and after removal are shown in Figure [12] It can be seen 
that after removal the problem of excess AM is eliminated 
and the profiles show a drop at low values of s, suggesting 
a > 1. These profiles are qualitatively similar to the ones ex- 
pected for exponential discs (see ^harma & Steinme tz ;2QQ5i; 



Ivan den Bosch et al.ll2QQlh . Although, high AM material is 
not removed but the curve at high AM end is still shifted to 
the left. This is because removing the low AM material in- 
creases the jtot of the system. In Figure [TS] we quantitatively 
assess the effect of removing the low AM material. We plot 
for various different values of sq, the mean fraction of mass 
removed /reject and the fraction of halos that can form expo- 
nential discs, i.e., halos with shape parameter a > 1.3. The 
higher the value of sq, the greater the amount of material re- 
moved and higher the fraction f{a > 1.3). The figure shows 
that if 90% of halos are to host exponential discs then one 
requires /reject to be 0.4. 
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Fig. 10. — The angular velocity ^2 as a function of radial distance for both 
gas and dark matter. The dashed line is the analytical expression that provides 
a good approximation to the profiles. 
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Fig. 12. — The effect of selectively removing low angular momentum ma- 
terial on the AMDs. The material was removed using Equation (|8) with 
5Q = 0.5. The solid curves show the mean AMDs computed over 42 ha- 
los and the dotted curves show the 16 and 84 percentile values. Note, jtot is 
not same for both the cases. In fact jtot increases after removal of low AM 
material. Also, for the after removal case, the area under the curve is less than 
one, the area actually represents the fraction of mass retained. 
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Fig. 11. — The shape in (Re, z) plane of the volume removed to get rid of 
the low AM material in galactic halos. The shape is controlled by the param- 
eter So which is the minimum angular momentum, as given by Equation ([8), 
that is retained. 

3.7. Spin Up and Spin Down of Halos Accompanied by 
Mergers 

It has been reported in earlier studies that immediately after 
the merger, e.g., the point of pericentric passage, the spin pa- 
rameter of the dark matter halo is found to be higher and later 
on as the system virializes the spin is found to drop. To study 
this we plot in Figure [14] the evolution of spin parameter of 
gas and dark matter for two merging scenarios, mass fraction 
fm = 0.5 (Sim-1) and /^ = 0.1 (Sim-7). Note, for other 
merging scenarios the results are quite similar to the Sim-1 
case. In both panels Adm is found to drop sharply from the 
point of first pericentric passage with a slow subsequent rise 
later on. The evolution of Acas is sensitive to the choice of 




Fig. 13. — The variation of fraction of lost mass, /reject, and the fraction 
of halos having a > 1.3 with the parameter so which controls the amount of 
material removed. 

fm but in general shows much less variation than that of Adm • 
The drop in Adm for /^ = 0.5 case is about 30% whereas for 
fm = 0.1 is about 80%. The sharpest fall of Adm is found to 
last for about 1 h~^ Gyr and occurs between t = 4 h~^ Gyr 
to t = 5 h~^ Gyr in the simulations. 

Next, we look at the evolution of the ratio jcas/JBU- For 
/m = 0.5, the gas has lost some AM by the time of the start of 
the merger, but after that the ratio jcas/JDM seems to remain 
constant. On the other hand, for fm = 0.1 case, the gas is 
found to gain AM from dark matter. 

Normally, the virialization ratio 2T//7 + 1 and the offset pa- 
rameter defined as Ar/rvir is used to detect such non relaxed 
halos. These are also plotted alongside. Here, T is the kinetic 
energy and U the potential energy of the system and the offset 
is defined by Ar = |xcm — x^ax den|- The most commonly 
used values of these quantities are —0.5 < 2T/U + 1 < 0.5 
and Ar/rvir < 0.1. It can be seen from Figure [14] that both 
these criteria have limited effect in detecting such cases. Our 
results suggest that a choice of Ar/rvir < 0.025 should be 
more effective in detecting such high spin systems. 

4. DISCUSSION AND CONCLUSIONS 
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Fig. 14. — Evolution of spin, total angular momentum, virial mass, offset 
parameter and the virialization parameter with time. Results are shown for 
mass ratio fm = 0.5 (left panels) and fm = 0.1 (right panels). A (10) and 
?7Zvir(10) refer to values of spin and virial mass at t = 10 Gyr h"-*^ and 
are used to normalize the values of spin and mass with respect to the final 
equilibrium distribution. 

We have performed non-radiative hydrodynamical simula- 
tions of mergers of spherical halos with a view to understand 
the AM properties of halos simulated in a cosmological con- 
text. The simulations being non-radiative are fundamentally 
aimed at determining the AMD of the gas before it cools onto 
the disk. In reality, at the time of most mergers a condensed 
component of gas also exists in addition to the hot gaseous 
halo. This cold component has not been taken into account in 
our simulations. This cold component plays a major role in 
fueling the growth of black holes and quasar activity. How- 
ever, a si gnificant fraction of t his is also driven away due to 
feedback (IHopkins et al.ll2QQ5ll2QQ6|) . The relative proportion 
of hot and cold component during a merger, their interaction 
with each other and their relative contribution to the final disc 
is still not fully understood. In this light, our results concern- 
ing properties of disc galaxies are mostly applicable for discs 
(or its parts) that are predominantly formed out of gas accreted 
from the hot gaseous halo. The main AM properties studied 
include the evolution of AM, the spatial distribution of AM 
and the orientation of the AM within the halo. We also ex- 
plored the differences between the AM properties of gas and 
dark matter and explain their origin. We now summarize our 
results and discuss their implications for the formation of disk 
galaxies. 

The shape parameter a of AMDs of gas in merger rem- 
nants is less than one for a wide variety of orbital parameters. 
This seems to be a generic result of the merging process. Val- 
ues greater than one and reaching upto 1.08 only occur for 
unrealistically large value of A or very low concentration pa- 
rameter. Lower values of mass ratio fm and higher values 
concentration parameter c result in lower value of a . Under 
the assumption that disks form under conservation of AM this 
leads to disks that are too centrally concentrated, as exponen- 
tial disks require a value of a greater than 1.3. In a previous 



study by iMaller & Dekell (120021) it was suggested that halos 
acquire most of their AM by means of major mergers while 
minor mergers with small satellites, which come in from ran- 
dom directions, contribute to the low AM material. They ar- 
gued that by preferentially discarding gas from the shallow 
potential of these small halos, e.g., by means of supernova 
feedback, the AMD problem could be solved. However, our 
results show that even in absence of minor mergers, the AMD 
generated by a major merger has an excess of low AM mate- 
rial. Even the most favorable of merging scenario thus cannot 
account for the formation of disk galaxies. Indeed, mergers 
in which the puffing up of gas by feedback was mimicked by 
decreasing the concentration parameter of the gas, did show a 
slight increase in value of a, suggesting it may partially help 
in reducing the low AM material but is not enough to solve 
the problem. 

We find that the angular velocity Q is almost independent 
of the spherical coordinate 6 but exhibits a significant ra- 
dial gradient. Such behavior is also seen for halos drawn 
from cosmological simulations. Hence, spherical shells of 
matter appear to be moving in solid body rotation. This 
seems like a safe assumption to be used for semi-analytic 
modelling (van den Bosch 2001, 2002). The spatial distri- 
bution of low AM material is found to be in the center and 
along a conical region along the rotation axis. This sug- 
gests that a mechanism which can preferentially eject the 
material in the center and along the poles can alleviate the 
AMD problem. In fact, feedback from intense star forma- 
tion in the inner regions can drive such an outflow . Addi- 
tionally, feedback from AGN jets is also expected to evac- 
uate material along the axis. Evidenc e of such conical out- 
flows is also provided by observations (Heckman et al." 1 9901: 
ShoDbell & Bland-Hawthorn 1998; Veilleux & Rupke 2Q02; 
Veilleux et all 120051: Bland-Hawthorn et al. 2007).^ More re- 
cently, using oxygen absorption lines in background quasars, 
significant amounts of diffuse metal rich gas has been de- 
tected in halos of star forming galaxies, suggesting large 
scale star formation dri ven galaxy outflows (iTumlinson et al.l 
1201 ll: PlVipp et al.ll201 1|) . Essentially, during the formation of 
the galaxy the star formation will be strongest in the cen- 
tral regions and this will drive an outflow which will ex- 
pand more along the rotation axis due to the flattened ge- 
ometry of the assembling gas. This will prevent the rest of 
the low AM material from fallin g onto th e disc. This mech- 
anism was initially discussed in ISharmal (|2005) (Sec-3 and 
5.2 of the PhD thesis) , and more recently has also been pro- 
posed by iBroo k et al I dToilb). Using high resolution cos- 
molo gical simula t ions inc orporating star formation and feed- 
back lBrooket al.l (|2011b|) clearly demonstrate that the above 
mechanism is resp onsible for the formation of bulgeless dwarf 
galaxies (see also 'Guedes et al.''201 1). Additionally, as sug- 
gested by Brook etal. (2011a), for small mass systems the 
outflows can eject the low AM material but for large mass 
systems it can also drive a fountain leading to mixing and re- 
distribution of the low AM material to be accreted later on as 
high AM material. In addition to cosmological simulations, 
preferential ejection of low AM material by feedback from 
the central regions was al so used in semi-analytic modelling 
by (iDutton & van den Bo sch 2009; Dutton 2009) to success- 
fully reproduce the exponential structure of discs. 

We tested a simple geometric criteria to selectively remove 
the low AM material from halos and found that resulting 
AMDs are in good agreement with those of exponential disc 
galaxies. Our results suggest that in order for 90% of ha- 
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los to form realistic exponential discs approximately 40% of 
the baryonic mass needs to be rejected. The presented cri- 
teria can be easily applied to dark matter simulations or in 
semi-analytical modelling. Note, our results provide an upper 
limit on the fraction of retained baryons fd = l — /reject = 
^disc/((^b/^m)^vir) = 0.6, Mdisc being mass of stars and 
gas in the galactic disc. Real galaxies can have even lower /d 
depending upon other physical processes which regulate the 
star formation efficiency. In the context of the missing baryon 
problem, the upper limit given by AMD argument is in good 
agreement with observations. iGuo et al] (i2010) using abun- 
dance of galaxies from SPS S find a maximum value of /a to 
be 0.2 at Mvir = 6 x 10^^ iMcGaugh et aH fcOlQ.) estimate 
fd < 0.4 for spiral g alaxies and even lower for dwarfs. Note, 
iMcGaugh et al.l (|2010) use an overdensity factor of A = 500 
(instead of 100) to measure the mass of the halo, so actual 
values of /d with in virial radius would be about a factor of 
two lower. 

The difference between collisionless dynamics and gas dy- 
namics results in differences between the AM properties of 
the gas and dark matter and this can potentially have implica- 
tions for studies that assume them to be same. The gas as 
compared to dark matter is more efficient in depositing its 
orbital AM in the central parts of the halo. This results in 
a higher value of the spin parameter A for the gas as com- 
pared to dark matter and also a moderately high value of the 
shape parameter a. Lower values of mass fraction fm and ini- 
tial concentration also result in higher Acas/ADM while lower 
values of Ainitiai and merging time torh result in lower values 
of spin ratio. About 6 h~^ Gyr after the merger, i.e., the first 
pericentric passage, the ratio Agas/AoM is found to be greater 
than 1 for all merging scenarios analyzed here. This seems to 
be consistent with spin ratios of halos obtained from cosmo- 
logical simulations at z = 0, where < Aga s/ApM > is close to 
1.4. Using a sample of 14 dwarf galaxies v an den B osch et al. 
{2001) had found the median spin of galaxies to be 0.06, as- 
suming Adm = 0.0367 this gives Agai/AoM = 1-63. The 
higher spin of galaxies could be due to the gas having higher 
spin, but it could also be due to preferential rejection of low 
AM material during the assembly of the disk. 

We find that for mergers with zero intrinsic spins, the AM 
vectors of gas and dark matter are well aligned with misalign- 
ment angle being less than 2°. On the other hand, mergers 
having non-zero intrinsic spins which are inclined at an angle 
to the orbital AM vector can result in a misalignment of about 
20°, consistent with halos simulated in a cosmological con- 
text. Since halos simulated in cosmological context undergo 
multiple mergers with different spin orientations, the above 
result provides a natural explanation for this. This shows that 
the misalignment can be explained purely by means of merg- 
ers without any need for the gas and the dark matter to be 
torqued differently during the formation of the proto-halo. In 
general, the gas within the virial radius is more effective in re- 
taining the information about the intrinsic spins of the merg- 
ing halos whereas the dark matter is more effective in retain- 
ing the orbital AM information. The misalignment between 
gas and dark matter has important implications for studies 
such as, the correlation between the anisotropic distribution 
of satellite galaxies and the major axis of the central galaxy 
and weak lensing studies attempting to measure the ellipticity 
of the dark matter halos. 

Mergers with non aligned spins also tend to make the AM 
of gas as measured in radial shells misaligned with each other. 
Since galaxies generally form inside-out, later infall of mis- 



aligned material can cause warps in disk gala xies, and this 
has recently been shown by lRoskar et al.l (1201 Q|) in cosmolog- 
ical hydrodynamical simulations with star formation. They 
find that immediately after the major merger the inner gas 
which forms the disk is misaligned with the rest of the gas 
in the halo. Later infall of misaligned gas causes the warps. A 
probable explanation for the cause of misalignment is given 
by them as the fact that interactions such as minor mergers 
can affect AM of the inner and outer regions differently. We 
have here explicitly demonstrated as to how major mergers, 
in which the orbital AM is not aligned with the intrinsic spin 
of the halos, generates such a misalignment. Minor mergers 
later on may further alter the orientation but are not necessar- 
ily required to generate the misalignment. 

Our results show that the orientation of the AM within the 
halo depends sensitively upon the orientation of the intrinsic 
spins of the merging halos with respect to that of the orbital 
AM. The larger the initial misalignment between the initial 
AM vectors the larger the final misalignment between the in- 
ner and outer parts. This suggests that warps may offer the 
possibility to probe the merger history of the halo. 

Observational eviden ce for w arps is quite ubiq- 
uitous dSancIsil [T976I: iBriggsl [l995 Eubi^ [1994]: 
IVerdes-Montenegro et all l2QQ2h . iGarcia-Ruiz et al.l (l2QQ2h 
find that in their sample of galaxies, all galaxies that have an 
extended HI disk with respect to the optical are warped. If 
misalignments in AM are not as frequent then this could pose 
a problem. Our results show Figure |8]that AM of gas within 
r < O.lrvir and r > 0.9rvir is misaligned by more than 10°, 
with respect to the total AM vector, for about 84% of the 
halos (the median misalignments being about 30° and 17° 
respectively). This demonstrates that the misalignments are 
quite common and supports the idea that they are responsible 
for warps. Just as perfect prograde and perfect retrograde 
mergers are rare so are systems with small angle warps and 
systems with counter rotating gas. In future, observations 
with detailed statistics on the orientation of the warps could 
be employed to check if they match with the distribution of 
misalignments predicted by theory. 

As mentioned earlier, the amount of misalignment in gen- 
eral is found to increase with the increase of angle between 
the orbital and intrinsic AM vectors. For retrograde encoun- 
ters the gas at intermediate radii is even found to be counter 
rotating. This could be responsible for the counter rotating gas 
seen in some galaxies. Generally, mergers of gas rich systems 
are invoked to explain such systems. The quantity of counter 
rotating gas in some galaxies such as NGC3626 is so large 
that a single minor merger cannot properly account for it. If 
on the other hand a merger is not minor then it can heat up and 
thicken the disk considerably. A slow, continuous and well 
dispersed accretion, as opposed to an accretion via a merging 
system is preferred (Ciri et al. 1995; Thakar & Ryden 19%|). 
Counter rotating gas in galactic halos formed by retrograde 
mergers as shown here, naturally provides such an extended 
reservoir of gas. A recent merger which can potentially heat 
up the disk is not required, the counter rotating gas is formed 
early on during the last major merger, which causes the inner 
and outer regions to rotate in different directions. In such a 
scenario, the inner regions first assemble to form the disc, rest 
of the material falls later on to generate the counter rotating 
gas. 

We also studied the issue of spin up of a halo 
undergoing a merger and t he su bsequent spin down 
during virialization (.Gardnerl I2QQ1I: .Vitvitska et al.l I2Q02I; 
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iPeiran i et alJl2004t iHetznec ker & Burkertll2Q06 h. As argued 
bylDfOnghia & Navarro (20Q7|), in coUisionless mergers cen- 
tral regions tend to be populated by low AM material and high 
AM material is pushed to weakly bound orbits. When AM is 
measured with the fixed radius, such as the virial radius, the 
effect is a spin down. Our merger simulations also show a 
similar effect. For dark matter the inner half loses AM while 
the outer half gains. The main cause for such a redistribution 
of AM is the coUisionless dynamics and is as follows. For 
dark matter during the collision the late in-falling particles 
have high AM and high energy and they gain energy during 
the collapse and hence can easily climb out of the final relaxed 
potential which is much shallower. Hence, high AM particles 
get pushed to weaker and weaker orbits making them move 
outwards. For gas the late in-falling particles shock and de- 
posit their AM onto the inner regions making them behave 
differently. We find that for a Milky Way sized halo the spin 
down process lasts about a giga year, and the spin of dark 
matter can fall by about 40 — 80% during this time, depend- 
ing mainly upon the mass ratio of the merging components. 
The spin of gas shows somewhat less variation. We find that 
the virial ratio 2T/U + 1 is not very effective in detecting 
such situations. The offset parameter is more successful in 
detecting such cases but a value of Ar/rvir < 0.025 would 
be needed, which is much less than what is currently used 
(iD'Onghia & Navarroll2007l: iNeto et al.ll2007h . Hence, recent 
results showing high spi n systems to be more clus tered may 
be aff ected by this bias (iBett et al.ll2"007l: iDavis & Natarajan 
|2010). Alternatively, it may reflect the fact that in clustered 
environments mergers and hence non relaxed halos are more 
common. If non relaxed halos is the cause of correlation 
between the clustering and spin then the ability to observa- 
tionally detect it by measuring spin of galaxies is unclear, as 



galaxies form out of baryons in relaxed halos. Additionally, 
it is not known if the spin of baryons would also show such a 
clustering. 

Finally, our results show that mergers of NFW halos nat- 
urally generate the universal form of AMDs as seen in sim- 
ulations. For dark matter the value of the shape parameter 
a is in excellent agreement with the results from cosmolog- 
ical simulations. However, this does not mean that merg- 
ers are the only way to generate such d istributions. As has 
been shown recently by 'Wang & Whitd (|2009) even hot dark 
matter simulations which have almost no mergers show such 
AMDs. Hence, the origin of the universal form is more gen- 
erally related to the virialization processes such as the violent 
relaxation. However, our results show that mergers do induce 
subtle differences between the AM properties of dark matter 
and gas. The alignment of the AM vector within the halo 
and also that of gas with respect to dark matter is sensitively 
related to the merger history and may serve to discriminate 
the dark matter models. In hot dark matter models, although 
less likely, misalignments as discussed above could also be 
produced if matter coming from different regions have AM 
pointing in different directions. 
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APPENDIX 
ORBITAL PARAMETERS 

The merger of two bodies of mass mi and m2 can be reduced to the motion of a test particle, with a reduced mass /i = 
mim2/{mi + 7712), in the potential of a mass M = mi + m2 (Figure [T5]). The initial conditions are set by specifying the relative 
separation Trei and relative velocity Vrei- In a cosmological context the two masses first move apart due to Hubble expansion 
and eventually, come to a halt and collapse due to their mutual gravitational attraction. The orbits of interest are those which are 
bound and collide within a Hubble time. A bound orbit can be fully characterized by its eccentricity e and the semi-major axis a. 
The energy of the orbit £^orb, and the orbital time period Torb are related to a by 



E, 



orb " 



GMfi 
2a 



, ^orb = 27r 



GM 



Eorh can be written in terms of Torb as 



Eor, = -l{4n'G'y/X~"J'UM'/' 



where f^ = ^- The angular momentum Lorb is related to eccentricity e by 



^ , /<;'/Vf5/2/-3/2 /-, _ 2 

V^l^orbl 



(Al) 



(A2) 



(A3) 



According to the tidal torque theory the system acquires AM during its ex pansion phase, with the AM incre asing nearly linearly 
with time during the initial linear phase of growth of density perturbations (Iwhitell 1 9841: iDoroshkevichll 1 97Q|) . The acquisition of 
AM ceases in the non-linear regime. We assume that all AM is acquired by the time of maximum expansion which gives 

rrei = a(l + e) (A4) 

At maximum expansion, the radial velocity being zero, the total velocity is given by the tangential velocity. 

Torb 



^rer 



/irrel 



(A5) 
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Fig. 15. — Merger of two halos can be reduced to a one body problem of mass /i moving in the potential of mass M. The orbit can be characterized by 
semi-major axis a and eccentricity e. At maximum separation r^^i = a(l + e) the tangential velocity v^^i is given by the angular momentum acquired by the 
masses during the expansion phase. 
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Fig. 16. — Dependence of fraction of mass lost during a collision on merger parameters. The fraction of mass lost is an increasing function of the kinetic 
eAergy involved in the collision and a decreasing function of the total binding energy .of the system. .. i 

Since the merging bodies are extended objects, the total energy is given by the sum of the orbit al energy plus the se l f ener gy of 

the bodies. The self energy of a body of mass M^ and radius R^ , having an NFW density profile (iNavarro et al.ll 19961 Il997h with 

concentration parameter c, is given byQ 



Ev — —fc 



GMl 



where fc 



2~ 



l/(l + c)2-21n(l + c)/(l + c) 



2Ry ''' 2 (ln(l + c)-c/(l + c))2 

Assuming that both halos are virialized at a redshift of z, My can be written in terms of Ry are as 



M.,=^A„)5|^ = «r 



2G 



(A6) 



(A7) 



where A{z) is the over-density criteria used to identify a virialized regi on, i.e., a spherical regio n whose average mass density is 
A(z) times the critical density at that redshift. A(z) is approximated by (iBryan & Normanll998l) A(z) ~ (187r^ + 82a;+— 39a;^), 
where x = ^m{z) — 1. Consequently, the total energy is given by i3 = Eyi + Ey2 + £'orb- Analogously, the spin parameter A 
of the whole system is given by 



A = 









1/2 



(A8) 



Instead of the semi-major axis a and the eccentricity e the orbit can be equivalently parameterized in terms of the orbital time 
period Torb and the spin parameter A. We restrict ourselves to values of Torb which have r^ei > ri2 where ri2 = r^iri + ^vir2. 



MASS STRUCTURE OF REMNANT HALOS 



vir of the remnant is less than 
mvir)/^tot- Also, the 



The final properties of the merger remnants are given in Table \T\ We note that the virial mass m 
the total mass of the system mtot- Hence, a fraction of mass is lost which we define as /lost = (^tot 
concentration parameter of the remnant halo Cfinai is slightly larger than cinitiai- 

It is interesting to know if the final properties of the halo e.g., mvir and Cfinai can be predicted from the initial conditions. 
We expect the fraction of lost mass /lost to be an increasing function of the kinetic energy KE involved in the collision and a 
decreasing function of the total binding energy of the system. We find that the following empirical formula, which satisfies the 
above conditions, fits the results obtained from simulations (Figure [16]). 



„ Maximum KE of collision at rgep = ^12/2 

/lost OC 
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For the Einasto 



|^tot| + |2Fl2| 
profile the formulas are available at 



where V12 = 
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INichols & Bland-HawthornI ( f200^.[20TTh 
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If Cinitiai is higher the system has higher |£^tot | consequently, it is more bound and loses less mass. If |£^orb | is higher the system 

is again more bound and also the KE of the collision is less, consequently reducing the mass loss. 

Interestingly, the total energy of the remnant halo £^vir (putting Cfinai and mvir from Table [T] in Equation (IA6I) ) is nearly equal 
to the energy of the system £^totai before the merger. This suggests that the mass that lies outside the virial radius, consists of a 
bound and an unbound part and has almost zero net energy. Consequently the concentration parameter of a remnant halo can be 
predicted from the knowledge of its orbital parameters. 
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